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Abstract. We study the behavior of nonzero rest mass spinning test particles 
moving along circular orbits in the Schwarzschild spacetime in the case in which 
the components of the spin tensor are allowed to vary along the orbit, generalizing 
some previous work. 
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1. Introduction 

The equations of motion for a spinning test particle in a given gravitational background 
were deduced by Mathisson and Papapetrou ^ |2] and read 

-i = - -R\ aP U»S a P ee , (1.1) 

dm 2 

— = P^U" - P"U>* , (1.2) 

aru 

where P M is the total four-momentum of the particle, is the antisymmetric spin 
tensor, and U is the 4-velocity of the timelike "center of mass line" used to make the 
multipole reduction. In order to have a closed set of equations, Eqs. (|1.1|) and l|1.2|l 
must be completed by adding supplementary conditions (SC) whose standard choices 
in the literature are the 

1. Corinaldesi-Papapetrou [3j conditions (CP): = 0, where the index v 
corresponds to a coordinate component and t is a timelike slicing coordinate, 

2. Pirani 4 conditions (P): S^ V U V = 0, 

3. Tulczyjew conditions (T): ' P v = 0. 
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Only solutions of the combined equations for which both U and P are timelike vectors 
are considered, in order to have a meaningful interpretation describing a spinning test 
particle with nonzero rest mass and physical momentum. 

Not much is known about actual solutions of these equations in explicit spacetimes 
which satisfy the Einstein equations. In a previous article [H], we considered the 
simplest special case of a spinning test particle moving uniformly along a circular 
orbit in the static spherically symmetric Schwarzschild spacetime, but because 
these equations are still complicated, we looked for solutions with constant frame 
components of the spin tensor in the natural symmetry adapted static frame, i.e., 
coinciding with a static tensor field along the path. Such a static spin tensor is a 
very strong restriction on the solutions of these equations of motion, leading to special 
solutions in which the spin vector is perpendicular to the plane of the orbit, and 
contributes to an adjustment in the acceleration of the orbit. 

Here we consider the slightly less restrictive case where the spin components are 
not constant, but the motion is still circular. However, in this case it is clear that 
if the spin tensor has time-dependent components, its feedback into the acceleration 
of the test particle path will break the static symmetry of that path unless the spin 
precession is very closely tied to the natural Frenet-Serret rotational properties of the 
path itself. Indeed we find that only the Pirani supplementary conditions permit such 
specialized solutions since they allow the spin tensor to be described completely by a 
spatial spin vector in the local rest space of the path itself. By locking spin vector 
precession to the Frenet-Serret rotational velocity of the path, solutions are found with 
a spin vector Fermi- Walker transported along an accelerated center of mass world line. 
The remaining choices for the supplementary conditions have no natural relationship 
to the Frenet-Serret properties of the particle path and do not admit such specialized 
solutions. 

With the assumption of circular motion, one can solve the equations of motion 
explicitly up to constants of integration. By a process of elimination, one can 
express them entirely in terms of the spin components and particle mass as a 
constant coefficient linear system of first and second order differential equations. 
By systematic solving and backsubstitution, one gets decoupled linear second order 
constant coefficient equations for certain spin components, which are easily solved 
to yield exponential or sinusoidal or quadratic solutions as functions of the proper 
time, from which the remaining variables may be calculated. Imposing the choice 
of supplementary conditions then puts constraints on the constants of integration or 
leads to inconsistencies. The details of the decoupling and solution of the equations 
of motion are left to the Appendix, leaving the imposition of the supplementary 
conditions to the main text. 

2. Circular orbits in the Schwarzschild spacetime 

Consider the case of the Schwarzschild spacetime, with the metric written in standard 
coordinates 



and introduce the usual orthonormal frame adapted to the static observers following 
the time lines 




(2.1) 



e t - - (1 - 2M/r)- 1/2 d t , e? = (1 - 2M/r) 1/2 d r , e § = -de, e 



r sin 



1 



9 



d<t>, 



(2.2) 
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with dual frame 

J = (1- 2M/r) 1/2 dt, J = (1- 2Af/r)" 1/2 dr , J = rd9 , ^ = rsin(9d</>, (2.3) 

where {9*, d r , dg, c^} and {di, dr, d#, d</>} are the coordinate basis and its dual, 
respectively. 

In order to investigate the simplest special solutions of the combined equations 
of motion, we explore the consequences of assuming that the test particle 4-velocity 
U corresponds to a timelike constant speed circular orbit confined to the equatorial 
plane 9 = tt/2. Then it must have the form 

U = T[d t + ^]= 1 [e i + ve^ 7 = (1 - v 2 )~ l l 2 , (2.4) 

where £ is the angular velocity with respect to infinity, v is the azimuthal velocity as 
seen by the static observers, 7 is the associated gamma factor, and T is a normalization 
factor which assures that U ■ U = — 1. These are related by 

C = (-9tt/g^) 1/2 v , r = (-9u - C 2 5^r V2 = {-9tt)- l,2 l , (2.5) 

so that £r = 7zV(fl'0<?O 1 / 2 , which reduces to ^vjr in the equatorial plane. 

Here £ and therefore v are assumed to be constant along the world line. We 
limit our analysis to the equatorial plane — 7r/2; as a convention, the physical 
(orthonormal) component along — dg which is perpendicular to the equatorial plane 
will be referred to as "along the positive z-axis" and will be indicated by the index z 
when convenient: = — es. Note both 9 = ir/2 and r = ro are constants along any 
given circular orbit, and that the azimuthal coordinate along the orbit depends on the 
coordinate time t or proper time r along that orbit according to 

4> ~ 4>o = (t = Qui~u , £lu = l v l r , (2-6) 

defining the corresponding coordinate and proper time orbital angular velocities £ and 
flu. These determine the rotation of the spherical frame with respect to a nonrotating 
frame at infinity. 

Among all circular orbits the timelike circular geodesies merit special attention, 
whether co-rotating (( + ) or counter-rotating (£_) with respect to increasing values of 
the azimuthal coordinate cj> (counter-clockwise motion). Their time coordinate angular 
velocities £± = ±(k = ±(M/r 3 ) 1 / 2 , which are identical with the Newtonian Keplerian 
values, lead to the expressions 



U± = 7ic[ e |± v K eA , v K 



M 



2M 



1/2 



IK 



2M 



r-3M 



1/2 



(2.7) 



where the timelike condition vk < 1 is satisfied if r > 3M. At r = 3M these circular 
geodesies go null. 

It is convenient to introduce the Lie relative curvature [7[ |H| of each orbit 

1 / 2M\ 1/2 Ck 



fc(iic) = ~9f lnV500 = — 1 = . (2-8) 

r \ r J v K 

and a Frenet-Serret intrinsic frame along U [Hj, defined by 

E = U , E x = e f , E 2 = 7 ^e t - + e^] , E 3 = e 2 (2.9) 

satisfying the following system of evolution equations along the constant radial 

acceleration orbit 

DU , x DEi DE 2 DE 3 

— = a(U) = kE, , — J- =kU + nE 2 , = —T\E\ , — -1 = , (2.10) 

ottj dT\j dT\j dT\j 
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where in this case 



2 / 2 

; 2r 2 2 1 7 \ v — U K) /■ 



Tl = ~ = _ k <V*tt" = "72— Ot ■ (2-11) 



1 dn 7 2 7 2 /y 

rTT^ - = _ "'(lie) — Y V = 2 ' 

27 dv y ' Y K Tk v k 
The projection of the spin tensor into the local rest space of the static observers 
defines the spin vector by spatial duality 

S fi = \vJ lS {e t ) a S lS , (2.12) 

where n a p-ys = ^—g^ap-yS is the unit volume 4-form constructed from the Levi-Civita 
alternating symbol e Q/ 3 7< 5 (e^g^ = 1), leading to the correspondence 

(S f , S 8 = -S* S*) = 6%, -5^, S^) . (2.13) 

For the CP supplementary conditions only these components of the spin tensor remain 
nonzero, while in the remaining cases the other nonzero components are determined 
from these through the corresponding orthogonality condition. The total spin scalar 
is also useful 

* 2 = = ~S* fi - Si - 5| t - + S% + + , (2.14) 

and in general is not constant along the trajectory of the spinning particle. In the 
Schwarzschild field the total spin must be small enough compared to the mass of 
the test particle and of the black hole \s\/(mM) <C 1 for the approximation of the 
Mathisson-Papapetrou model to be valid. This inequality follows from requiring that 
the characteristic length scale |s|/m associated with the particle's internal structure 
be small compared to the natural length scale M associated with the background field 
in order that the particle backreaction can be neglected, i.e., that the description of a 
test particle on a background field make sense |1(J|. 



3. Solving the equations of motion: preliminary steps 

Consider first the evolution equation for the spin tensor l|1.2|) . By contracting both 
sides of Eq. I|1.2|) with U V) one obtains the following expression for the total 4- 
momentum 

DS^ 

P" = -([/ • P)U^ -U v — ^mW+P?, (3.1) 

dry 

which then defines the particle's mass to, which a priori does not have to be constant 
along the orbit, while Pjf = U a DS afi /dry is the part of the 4-momentum orthogonal 
to U. Finally let U p denote the timelike unit vector associated with the total 4- 
momentum P = \ \P\\U P = fi U p . 

Backsubstituting this representation Eq. (|3.1|) of the momentum into the spin 
evolution Eq. (|1.2|) expressed in the static observer frame leads to 

n d ^ dS*,;,-. Ck , 2 2 \ c /o n\ 

0= d^r-^ +7 ^ ( 3 - 2 ) 



_ ^ jjjg _ 7^ Of 
drrr drr " 1 



v-P^-^-Ste, (3.3) 



W arjj 7^. v K 
^-7^Ck^ + 7 ^%. (3.4) 
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From i|3.1[l . using the definition of P s and equations i|3.2[l - i|3.4[l . it follows that 
the total 4-momentum P can be written in the form 



P = 7 (to + vm s )ei H 



7 



dru 



Ck q 

■ 7^ — <% 

V K v 



1 

Bf H — 
7 



+ 7(mz/ + TO s )e^ 
to?7 + m.-Bi H — 



7 



dS^ 
drr/ 



Ck a 



t<f> 



Gf H 

7 



dS'ffl 

7^-7^ 



dri 



r./ 



with 



(3.5) 



(3.6) 



dry zyj^ 

Next consider the equation of motion . The Riemann tensor spin-curvature- 
coupling force term is 



^ (sc) =7C^%e t - 



2S if 



s ie + 2vS §i 



S t4> e i\ 



(3.7) 



Using H3.1fl . the balance condition which allows a circular orbit of this type to exist 
can be written as 

ma(U) = F {so) + F {sc) , (3.8) 

where a(U) is the acceleration of the [/-orbit and F^ so ^ = —DP S / ' drrjj defines the 
spin-orbit coupling force term, which arises from the variation of the spin along the 
orbit. 

Taking i|3.5fl and i|3.6|) into account, Eq. (|l.lf> gives rise to the following set of 
ordinary differential equations 







d 2 S t 



t<ji 



dr 2 



I^vvkCk- 



dS A 



dru vk 
l 2 vC 2 K (v 2 ~v%)S id T - 
d 2 S„. < 12 



7— +v K )- 



r<p 



C: 



d 2 ^ 

dr 2 
dr* 



7 2 ^ 2 n 
ITk v k 
d 2 S &1 



t( t> ' dr v 

Ck AS ij> 
7^ v K dru 



dru 



+ vYCk(v 2 -v k )S^ 



S, 



dr 2 



Ck , 2 
7 — i/ 

VK 



VK 

4)^ + 2Ki^ 



(if / 2 2i dS 1 ^ 
7— (f + z^)- 



VK 



dru 



2jv 



dru 
Ck dSf f 



CK^fa , 



vk dru 



Ck 



2 2 

Tk v k 



dm 
dr v 



(3.9) 

(3.10) 
(3.11) 

(3.12) 



Note that there are two equations containing the second derivative of SW; this is due 
to the presence of its first derivative in two different components of P (more precisely, 
in P* and P$, see Eqs. (jHSl and fity ). 

Once the system of constant coefficient linear differential equations H3.2jl - I|3.4|l 
and H3.9|) ~ H3.12|1 is solved for to and the spin tensor components, one may then 
calculate P. The system must be decoupled, leading to functions which are either 
exponentials, sinusoidals, or at most quadratic functions of the proper time along the 
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particle world line. The elimination method for decoupling the equations is crucially 
different depending on whether v has the values or ±vk or none of these values, 
since one or the other or neither term drops out of the spin equations (|3.2|) and so 
must be considered separately. From the details of their derivations discussed in the 
Appendix, one sees why there are several zones approaching the horizon where the 
solutions change character. 

4. Particles at rest: the v = case 

When the particle is at rest, the solutions for the components of the spin tensor and 
the varying mass m of the spinning particle are given by 

(i) 2M < r < 3M: 



(ii) r = Ml: 



(iii) r > 3M : 



S-a 



c 2 e 



c 3 e 



VK 



2c r) 



2+4 



= c 4 e wnr + c 5 e^° T 

JjJOT I „ .-WOT 



Sfj, — 



c 6 e~"- +c 7 e 
IkVk [c±e aoT - c 5 e-° 0T ] + c 8 
Ikvk \c^ ar - c 7 e-° 0T ] + c 9 





= Cl , 














= c 2 e7/< 3M ) + 


c 3 e 


r/(3M) 


+ V3Mc m 




m 


A T 


/(3M) 


4- c 3 e" 


-r/(3M)j + 


2 




= c 4 r 4 


■c 5 , 












= c 6 r 4 














_ V3 
9M 


C 4 2 

.y r " 


f c 5 r 


4- c 8 








_ V3 
9M 


"C6 2 

.~2~ T " 


f C7T 


4- Cg 







(4.1) 



(4.2) 



Cl 



% = c 2 e^ + c 3 e-^ r 



Ck2 + u 2 k ' 



2c n 



- VKtK[c2e" lT + c 3 e~^ T } , , . 



A" 



^fe = c 4 cos wqt 4- C5 sinwor 
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^1$ = c 6 COSUJ o T + c 7 sin loqt , 
S f g = IrVk [C4 sin uj t - c 5 cos uj t} 
S-j, = Jkvk [ce sinw T - C7 cosw t] 
where c m , ci, . . . , eg are integration constants and 



c 8 
c 9 



w = iujo = 



f M(2r-3M) 
r 3 (r-2M) 



From Eq. I|3.5[) the total 4-momentum P then has the value 



P = me { + Ui[c 2 e UlT - 
in cases (i) and (hi), and 



let 





C8 






a c 9 




6 ll- 




VK 


* 7|. 



p 



""•'•d\7- r; 



r/(3M) . 



r/(3M) 



V3 Q 
9M f « 



c 4 



V3 
9M 



5 f i - c 6 



(4.3) 



(4.4) 



(4.5) 



H ( 4 -6) 



in case (ii). 

At this point the supplementary conditions impose constraints on the constants 
of integration which appear in the solution. For a particle at rest (y = 0), the CP and 
P conditions coincide and imply that = 0, namely 

C2 = c 3 = c 4 = c 5 = c e = c 7 = , c m = , (4.7) 

leaving arbitrary values for ci , Cg, eg. As a consequence, m should be as well, implying 
that P should be spacelike and therefore physically inconsistent. 
The T supplementary conditions when v = imply instead 



i§ dr v 



dS: 



Stx 



to 



t<t> 



dru 
dS ie 



if dr v ' 
dSu 
^ dru 

q dS if 

g _ dSff 



q ^jt£ 

^ d Tu 



v k C,k[S^S^ + S i§ S~ § ] 



VK(K{S.g + S^)-mS if , 



vkCkS^S^ - mS i§ , 
vkCkS-qS^ + mS { £ . 



(4.8) 



- - dr v r< P dru 

By substituting the solutions given by Eqs. (|4.1|) - (|4.3|) into these equations, one finds 
that all the integration constants except c\ must vanish. This in turn implies m = 0, 
which again leads to a spacelike P. Thus a spinning particle with nonzero rest mass 
cannot remain at rest in the given gravitational held. 



5. Geodesic motion: the v = ±vk case 

When the test particle's center of mass moves along a geodesic (the orbit has zero 
acceleration a(U) = 0) with azimuthal velocity v — ±vk, the spin-curvature and 
the spin-orbit forces balance each other (see Eq. EHJ): F (so) = -F( sc ). The 
solution of Eqs. I|3.9[l - i|3.12[) determines the spin which leads to this balancing. In 
the Schwarzschild spacetime, timelike circular geodesies only exist for r > 3M. We 
consider separately the various cases: 
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S§4> = c 3 cos W 4 T + c 4 s i n W 4 T ) 

Oia = C5 COSW3T + C6 sinu^r ± — 11 , 

S f § = ^kvk [C5 sinw 3 T - c 6 cosw 3 r] , 
^ — c m , 



& = c 7 e° 2T + c 8 er^ T ± 2 



7k c 2 
Ck4-3 7 | ' 



f0 



± z^[c 7 e u 



c 8 e 



Ck '4 - 3 7 



2" ' 



S# = ± — (3 7 | - 4)- 1 /2[ Cse -^- _ C7e ^-] + C9 _ 3 7 | 



C2 



if 



(ii) r = 6M: 



i = C3 cos 7777-7 + C4 sin ■ 



18M 



18M 



\/6t . \/6i 



S » =C5COS 3M7 + C6Sill 36M ±2 ^ 



^ - 3 



C5 sin 



36M 



C6 COS 



V6r 
36M 



771 C m , 

Sf f = C 7 T + C 8 , 



% = T 



V2r 
12M 



C7 



"T + C 8 



C2 ; 



(iii) r > 6M : 



= c 3 cos 0J4.T + c 4 sin ui 4 t , 
Six = C5 cos UJ3T + cq sinter ± — - 



vk 



S f § = Ikvk [05 sinw 3 T - c 6 cosw 3 r] , 
nfi — c m , 



«SV f = c 7 cos cj 2 t + c 8 sina> 2 r ± 2 



7 K c 2 
4 - 3 7 | 



c _i_ r , 1 I , r, V K c 2 

= ± vk [C7 cos w 2 t + c 8 smw 2 Tj + c x + 2— 7ff - 



Ck '"4-37 



2 ' 
2 

&i = ± — (4- 3 7 | : ) _1/2 [c 8 cosa; 2 T- C7sina) 2 T] +c 9 - 3-7^ 

v IK 



c 2 



4-3 7 
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where c m , c\, . . . , Cg are integration constants, and three real frequencies are defined 
for each open interval of radial values by 



w 2 = iu 2 = Of (4 - 3-f K ) 



2 xl/2 _ 



IM(r-6M) 



r 3 (r-3M) 



I M 

^3 = (k = -jT 



1/2 



lo 4 = iu>i = Cx(37i - 2) 1/2 = -J- 



M 



,-3M- M 

Consider first the open interval cases r 6M. From Eq. I|3.5[l . the total 4- 
momentum P is given by 



P = <j m^K T C-ff 



S fi - (1 - 2^) 7 |ci - 2. 

(-2 



VK 



C2 



(1 + 2*4) 



t<t> i _ 4^2 



+ (1 - 4i^)c9 !> e f 



K 



\ —S-n T (1 + 2^|-) 1 ^ 2 [03 sinw4T — c 4 coso^t] 1 ■ 



i _ Of 

± \ mj K v K =F 



S f i - (1 ~ 2^) 7 |c! - 2 



vk c 2 
7x 1 - 4v| 



(5.5) 



We next impose the standard supplementary conditions. The CP conditions imply 
that Sfa = 0, namely 

c 2 = c 3 = c 4 = c 5 = c 6 = c 7 = c 8 = c 9 = , (5.6) 

so that the only nonvanishing component of the spin tensor is S z = = c\ = s, 
leaving arbitrary values for c TO as well. From Eq. (|5.5|) . the total 4-momentum P 
becomes (using m s = sjvk(k which follows from Eq. fi.ty ) 

P = mP± + sjvkCkE^ , (5.7) 

with U± given by Eq. I|2.7[) . Re-examing Eq. (|3.7|l shows that the spin-curvature force 
then acts radially, balancing the radial spin-orbit force. 
The P conditions imply 



% = > 



5 rf ± vkS+i = 



which lead only to the trivial solution 



Sn f ± vkS&i = , 



Ci = c 2 = c 3 = c 4 



C6 = C7 = C 8 = Cg = , 



(5.8) 



(5.9) 



with c m arbitrary, or in other words the components of the spin tensor must all be 
zero, which means that a non-zero spin is incompatible with geodesic motion for a 
spinning particle. 

The T supplementary conditions when v = ±i>k imply 



= V,..^ + S f A + -fcsJp*- ± tk&vkS. 



J te 



t<f> 



o = &,- d5 » 



il[±v K s if -s fi Y Sh ' 



Of r a 2 2 o2 1 

■^K — [S t .-v K S f§ \ 



dru '" l tr r <t" dr v 

S% T v K {\ + v 2 K )S if S^ + v'kS'^I - rni 2 K [S if T v K S f ^ , 



^3 Of 
Ik — 
vk 



vk 

.2 c2 
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= 7if[±^%-%] 



0= : - i«,ssJ S 



dS t$ 

dm 



■S 



;•(> 



dn, 



Ck „ 
lK—t>i f 

VK 



c 



J f8 u t4>J ' 



dm 



q d% 
fl * dru 



q ^% 



VK v v 



t<t> 



(5.10) 



By substituting into these equations the solutions given by Eqs. (15. 1|1 — (|5.3|1 . we obtain 
the conditions 



(5.11) 
(5.12) 
(5.13) 



C2 = c 3 = c 4 = c 5 = c 6 = c 7 = c 8 = , 
implying that the only nonvanishing components of the spin tensor are 
S z = -S P £ = -a , = c 9 , 

and either 

ci,c 9 arbitrary, c m = ±j k (kCi , 

which implies that spin component S z is proportional to the mass, locking them 
together by a constant of proportionality depending on the orbit velocity, or 

c\ = = eg , c m arbitrary , (5-14) 

corresponding to the zero spin case where geodesic motion is of course allowed. In the 
former case the total spin invariant (|2.14l) reduces to 

s 2 = -4 + c\ , (5.15) 

so that the condition \s\/(mM) <C 1 preserving the validity of the Mathisson- 
Papapetrou model reads 



1 



1- ^ 



1/2 



< 1 , 



(5.16) 



implying either c\ > eg or r ~ 3Af (where — > 00). In the limit r — > 3M where the 
circular geodesies become null and require a separate treatment, one has a solution 
for which this spin component S z is fixed to have a value determined by the constant 
mass m and the azimuthal velocity, the t-<j> component of the spin is arbitrary. If 
one takes the limit m — > 0, then the component of the spin vector out of the orbit 
vanishes, leaving the spin vector locked to the direction of motion as found by jll) 
who discussed the null geodesic case using the P supplementary conditions, the latter 
being the only physically relevant in such limit. 

Finally consider the remaining case r = 6M. Eq. i|3.5|) then shows that the total 
4-momentum P is given by 



I 18M fe 6M 



± 5 




c 7 > e? 



V6 
108M 



(5.17) 
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Imposing the standard supplementary conditions gives rise to the same result as for 
the general case r 7^ 6M. 

The CP conditions imply 

C2 = C3 = C4 = C5 = C6 = C7 = c 8 = , (5.18) 

so that the only nonvanishing component of the spin tensor is S z — —S f i — —c\, for 
arbitrary values of c m , leading to constant mass m. The P conditions give only the 
trivial solution 

c\ = c 2 = c 3 = c 4 = c 5 = c 6 = c 7 = c 8 = , (5.19) 
with c m arbitrary leading to constant mass m. Finally the T conditions imply 

C3 = c 4 = c 5 = c 6 = c 7 = c 8 = , (5.20) 
so that the only nonvanishing components of the spin tensor are 

S* = S f i = -a , S i$ = c 2 , (5.21) 



and either 



or 



V2 

ci,c 2 arbitrary, c m = ± J^ c i (5.22) 



ci = = C2 , c m arbitrary , (5.23) 

with constant mass m in both cases. In the former case the spin invariant (|2.14|l 
reduces to 

s 2 = -4 + c\ , (5.24) 
so that the condition \s\/(mM) <C 1 reads 

18 / c 2 \ 1/2 

«1 , (5.25) 



mM y/2 
implying ci > c 2 . 

Thus if the center of mass of the test particle is constrained to be a circular 
geodesic, either the spin is forced to be zero or have an arbitrary constant value of the 
single nonzero component S z of the spin vector out of the plane of the orbit. 

6. The general case: v 7^ and v 7^ ±vk 

For general circular orbits excluding the previous cases v = and v = ±vk, the 
solutions of the equations of motion for the components of the spin tensor and the 
mass m of the spinning particle are 

Sg^ = A cos VLt + B sin tor , (6.1) 
Sv-CcceKT + DvntHT + FSH, F = ^ + 2 ^\ {1 _ ^ yy (6-2) 



7 Kk(1 + 2^)(^ 2 -^) 

nv K [v*{\ + 2vi) -4(1 

— vk!k[D cos£l\T — Csinfiir] , (6.3) 

2 2 
v K v l - v z K 



v K y 2 {l-^ 2 K ) + v 2 K {2 + v 2 K ) 



> j/2 ., . 2 , lVC m -~l K - C ° 

Of % + *4(2 + *4) L Ck \v-Vk) 
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C\e 



Sl+T 



c 2 e 



-n 4 



c 3 e 



Sir 



m 



vk 



Ck ^+4(2 + 4) 

■{(34 + $) [ Cl e n + T 
1 ^ 1 



c^e 

2 



(6.4) 



-7(1/ - v K )c m + 



v—c 
Of 



1 



2i/4(l + 24) 



c 2 e-°+ T ] + (34 - $) [c 3 e°- r + c 4 e 



-n_ri 



3,| 



7 0? v 2 

+ n 



2^(1 



24) 



[ci< 



c 2 e 



-n 4 



24(1 + 24) 



[c 3 > 



C4e 



-n_ri 



}, (6.5) 



(6.6) 



(6.7) 



where A, B, C, D, c m , cq, ■ ■ ■ , C4 are integration constants, and the real positive 
frequencies f2 and Oi are given by 



n = 1 Q K {l + 2 VK ) 1 ' 2 ^ , 1 = 



(6.8) 



assumed to be distinct for the above equations to be valid, and the remaining 
abbreviations are 

1/2 



Ik v K 



v 2 ±ik§ 



K 



l -~ 2 



1/2 



with 



2 9 
Tk v k 



(1 + 24) 



9 2 



(6.9) 



(6.10) 



8 (1 + 24) • 

The behaviors of the azimuthal velocities &, v and vk as functions of the radial 
parameter r/M are compared in Fig. ^ They all coincide at r = 6M, where 
v = v = Vk = 1/2; for 2M < r < 6M it is i> < D, while v > D for r > 6Af. 

The quantities f2± also lead to angular velocities for certain intervals of values 
of the azimuthal velocity v. In fact we are interested in those values for which Q + 
and/or f2_ are purely imaginary, since the imaginary parts can be interpreted as 
additional frequencies characterizing spin precession. One must distinguish the cases 
2M < r < 6M and r > 6M, referring to Fig. [T] and to Eq. l(0)l : 

a) r > 6M: 

- if v > v (region I), the quantities Q± are both complex; 

- if 1/ = 0+ = Q_ is purely imaginary, since v > D; 

- \i v < v <v (region II), Q _ is purely imaginary, while f2 + can be either real 
(even zero) or purely imaginary; 

- if v < v (region III), Q + is purely imaginary, while can be either real 
(even zero) or purely imaginary; 

b) 2M < r < 6M: 

- if v > v (region IV), the quantities f2± are both complex; 

- if v = v, = f2_ is real, since v < D; 

- if v < v (region V), f2+ is real, since v < v, while f2_ can be either real (even 
zero) or purely imaginary. 
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UN, 




\\ w 

\\ 




V \\ 






V i 








in 



10 12 



r/M 



Figure 1. The azimuthal velocities u (y^ in the plot), i> (u t in the plot) and vj^ 
(dashed curve) are plotted as functions of the radial parameter r/M. The dashed 
vertical line corresponds to the value r/M = 6, where they all coincide. There are 
five different regions (explicitly indicated in the plot), depending on the relative 
sign of the azimuthal velocity v with respect to v and u, which correspond to 
intervals for which Q± are real, purely imaginary or neither, as explained in the 
text. 



All of these remarks so far assume that the two frequencies and Oi are distinct, 
necessary for the decoupling procedure which leads to this solution. A different result 
follows in the special case = Sl\. This occurs for the particular value of the azimuthal 
velocity 



! M(r-3M) 



lK y K> V r(r-2M) ' 



(6.11) 



which vanishes at r — 3M and is real for r > 3M, rising to its peak speed at 
r ~ 3.934M and decreasing asymptotically towards the geodesic speed from below 
as r —* oo. The solutions for the components Sg^, and S-g of the spin tensor are 
given by 

Sxi = A cos fir + B sin Qt , 



J f9 



C-- 



3 70^0^2 



2 n 



cos Q,T 



D-- 



3 7o^o 
2 ft 2 



Ck At 



sin fir , 



n 2 



ZWkCAB + A ^t) + 2-2 (u B - v 2 K D) 

n 2 l 

3is v%(k(2A + BQt) - 2 — {v Q A - v\C) sinQr 

7o 



cos Or 



with 



ft = fii = 



2v% 



1 + + 



1/2 



3M 



r(r - 3M) + 3M 2 



1/2 



(6.12) 
(6.13) 

(6.14) 
(6.15) 



Spin precession in the Schwarzschild spacetime: circular orbits 



14 



while that corresponding to the remaining components as well as to the varying mass 
m are obtained simply by evaluating the general solutions (|6.4|) - (|6.7|) at v = u a . 
The reality properties of the quantities Q± can be determined as done in the general 
case, noting that uq < v (corresponding to region V) always holds in the interval 
3M < r < 6M. For r > 6M however, we must distinguish two different regions, a) 
6M < r < ro, with f = 6M(1 + y/2/2) w 10.24M such that u = v, where v < v 
(corresponding to region III), and b) r > ro, where vo > v (corresponding to region 
II). 

The behavior of S, U and P along the world line itself is completely determined 
by the initial conditions 



dS.A 



(6.16) 

T = 



drjj 

and the corresponding conditions on the mass m of the particle which follow from 
Eq. H6.7fl . Thus in the special case in which the "center of mass line" is directed 
along a circular orbit, the completion of the scheme for the spinning test particle is 
equivalent to a choice of initial conditions. 

In principle the components of the spin tensor which are not constants should 
precess with the different frequencies which appear in Eqs. H6.1(l - (|6.6(l . leading to 
non-periodic motion, a feature that seems to characterize the general situation in the 
Schwarzschild ^5] and Kerr |131 1141 H5) spacetimes. However, this does not occur 
in practice once the CP, P and T supplementary conditions are imposed, as we will 
see below. It turns out that the nonvanishing components of the spin tensor are all 
constant in the CP and T cases, while the motion is periodic with a unique frequency 
in the P case. As one might expect, the particle mass m turns out to be constant in 
all three cases. 

6.1. The CP supplementary conditions 
The CP supplementary conditions require 

S tf = 0, % = 0, % = 0. (6.17) 
From Eq. (|6.5|l this forces 

ci = c 2 = c 3 = c 4 = , c = - — {v 2 - v 2 K )c m . (6.18) 

V V K 

Substituting these values into Eq. (|6.4|) then gives 

Cm 

so from Eq. (|6.7|l we get 



Cm = 2 S r4, > 6 ' 19 ) 



S H = a=—-±-. (6.20) 
Finally, from Eqs. (|6.1() and (|6.3() it follows that 

% = - S r§ = • (6-21) 

However, from Eq. (|3.5|l it follows that 

P = -8V K <; K e } , (6.22) 

since m s = —sjvkCk = —m/v, a consequence of Eqs. (|3.6|l and l|6.20|l . This result is 
unphysical since the total 4-momentum P is spacelike. 



Spin precession in the Schwarzschild spacetime: circular orbits 15 

6.2. The P supplementary conditions 

The P supplementary conditions require 

S^ = 0, S fi + S^ = 0, S§t + S § ^ = 0. (6.23) 

Under these conditions the components of the spin vector Sjj in the local rest space 
of the particle, — ^i] a l3 ~ l5 U a S- f s, expressed in the Frenet-Serret frame, are just 

{Sh,S%,Sb) = i^SgfrS.g^S^) . (6.24) 

Comparing the hrst Eq. (|6.23|1 with Eq. (|6.6(l we get 

ci = c 2 = c 3 = c 4 = , (6.25) 

so that Sfo, S-^ and the particle mass m are all constant. Eqs. Ij6.4|l and 1)6. 5|l together 
with the second of the Pirani conditions Eq. I|6.23[l imply 



vv K 7 | ;/ 2 (2-5^) + 4(l + 2^) 



c 0=~~ ~75- T^To E75 ^ I !i n I o..2 x c ™ > ( 6 - 26 ) 

hence from Eq. I|6.7|) 

1 ^ 2 - ^ 2 

Cm = i 1 + ^ 2 (l-4*4) + 4(2 + 4)J m ' (6 ' 2?) 
Next by substituting these values of the constants cq and c rn into Eq. (|6.4|l . we obtain 

isfj = S* = S H = -^Yf^ri—^—r ■ ( 6 - 28 ) 

Finally comparing the last of the Pirani conditions Eq. Q6.23I1 with Eqs. 16.lfl - H6.2fl 
and Eqs. (|6.12() - (|6.13(l leads to two possibilities: either 

case PI: A = B = C = D = , (6.29) 

which places no constraint on v and the spin vector is constant and out of the plane 
of the orbit, or 

case P2: C = = D , F = v , (6.30) 
the latter of which (again from Eq. (|6.2|l ') leads to the special azimuthal velocity 
!- *4/ 4 \ V2 {M{r-9M/4)\ 1/2 



V = V ^= ±2VK \-lTi^) = ±2 { r\r-2M) ) ' ^ 

The case PI has been already considered previously 0, leaving only P2 to be 
considered here. The corresponding azimuthal speed u/p2) vanishes at r = 9 /AM 
and is real for r > 9 /AM, rising to a maximum speed of 1 at r = 3M, corresponding 
to the two null circular geodesies, and decreasing asymptotically towards twice the 
geodesic speed from below as r — ► oo. 

The corresponding values of 7 and f2 are respectively 

(1 + 2*4) 1/2 yMr-2M) 
^ 2 ) = - r -3M ' (6 ' 32) 

and 



= 7lCx[(4 - 4)(1 + 2.I)] 1 / 2 = ^ V4r ^ , (6.33) 
v ' r r — SM 
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using Eq. I|2.6[l . To get the anglular velocity of precession with respect to a frame which 
is nonrotating with respect to infinity, one must subtract away the precession angular 
velocity = jv/r of the spherical frame. In the case P2 one finds f2(p2) — = 0, 
so the spin does not precess with respect to a frame which is nonrotating at infinity. 
Substituting these values back into Eq. (|6.28() then leads to 



<?3 =± ^_. (6.34) 



The remaining nonzero spin components (|6.1|H(|6.3|) can then be expressed in the form 
S u = S f§ = 7 _1 [^ sin ^(P2)T!7 - Bcosn^Tjj] , (6.35) 



sh(o) = 

leading to 



S' 2 



(0) 



Sl(0) = ± 



2Q 



(P2) 





/ COS 


sin0 






1)- 


— sin 


cos (j) 


9 




V o 










The spin invariant (|2.14() becomes in this case 



1 



A 2 + B 2 



1 



4Ci4- 



(6.36) 



(6.37) 



(6.38) 



5if if - 

The Mathisson-Papapetrou model is valid if the condition \s\/(mM) <^C 1 is satisfied. 
From the previous equation we have that either 7 — > 00 or the sum of the bracketed 
terms must be small, i.e., 



A 2 



l M 2 



< 1 



B 2 



m 2 M 2 



< 1 



[4A/ 2 C|(4-4)]- 



< 1 . 



(6.39) 



The latter possibility cannot occur for any allowed values of r/M, since the third term 
(which is dimensionless) of (|6.39(l is always greater than « 1.88, as is easily verified. 
The former possibility is realized only in the case of ultrarelativistic motion, which 
Eq. (|6.32|) implies occurs only as r — > 3M, where the orbits approach null geodesies 
and the limit m — > forces the component of the spin vector out of the plane of the 
orbit to vanish, locking the spin vector to the direction of motion exactly as discussed 
by Mashhoon [TTj . 

It is well known that the spin vector Sjj = S\jEi lying in the local rest space of 
U is Fermi- Walker transported along U in the P case, so it must satisfy 







D(iw)Su 
drr/ 



P(U) 



drr/ 



dr f , 



s u T i 



E 1 



dSl 
dru 



s b T i 



E-> 



(6.40) 



from (|2.1U|) . where P(U)£ = 5£ + U^U a projects into the local rest space of U. To 
check this we must show that the following two equations are identically satisfied 



dSh 







dS 2 



nSh = 



(6.41) 



dru ' dru 

But these two equations follow immediately from (|6.35(l . since t\ — f2(P2) results from 
the direct evaluation of the expression (|2.11(l for n, with v given by Ij6.31|l . 
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Thus given the rest mass m of the test particle, the constant component of the spin 
orthogonal to the orbit is then fixed by the orbit parameters, while the component 
in the plane of the orbit as seen within the local rest space of the particle itself is 
locked to a direction which is fixed with respect to the distant observers, since the 
angle of precession with respect to the spherical axes is exactly the azimuthal angle 
of the orbit, but in the opposite sense. In other words the precession of the spin, 
which introduces a time varying force into the mix, must be locked to the first torsion 
of the orbit itself in order to maintain the alignment of the 4- velocity with a static 
direction in the spacetime, and the spin does not precess with respect to observers 
at spatial infinity. Furthermore, the specific spin of the test particle cannot be made 
arbitrarily small except near the limiting radius where the 4-velocity of this solution 
goes null, and the spin vector is then locked to the direction of motion. Apparently 
the imposition of a circular orbit on the center of mass world line of the test particle 
is just too strong a condition to describe an interesting spin precession. 

The total 4-momentum P given by Eqs. I|3.5[) and (|3.6[1 can be written in this 
case as 



P = mU + m s E 2 + 7~ 



dru 



= mU + m s Ei — 



7^ 
r 



vkC,k Sjjez 



with v given by (|6.31|) and m s a constant 

Of ,_.c ..2 o A 



l—( vS tr 



'K 



2 Of , 2 

7 — \v 

VK 



v 2 k)S { 



(6.42) 



(6.43) 



but the final term in P (out of the plane of the orbit) oscillates as the spin precesses 
in the plane of the orbit. Note that the radial component of P is zero. 
The spin-curvature force (|3.7[1 simplifies to 



7&{ 



2S f ^ + vS, 



f(f) 

37 2 i<|(Sf}e f - S^eg) 



2vSz 



while the term on the left hand side of Eq. Ijl.lfl can be written as 



DP 

= [mK — m s Ti\ef 

dru 



S u +2vS, 



The force balance equation (|3.8|) reduces to 
ma(U)r = F^ so) + if c) , 
ma(U) d = = Fi 



(so) 



Fi 



(sc) 



(6.44) 



(6.45) 



(6.46) 



where 



l(U)r = 



F. 



(so) _ 



p( sc ) _ 



<p 

dru 

(so) 



= rn s Ti 



(6.47) 
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n <? dS te , c dS if , 2 o dS ^ 



u 



+ mj 2 vS. 



t<t> 



7^{h%<% + VKStgSf}] - j 2 S { ^[vS if - !^%]} 



= &, a 



r9 dr. 



u 



dr f , 



.Or 



VK 



4 - „(i + v\)S if S fi + vlS), - mi 2 [S if - uS H ] 



i Ck a 



+ ^d^7 



_ 7 dTt; + dr v + b * d W + 7 v K V 6 * ™ K ^ 1 



(6.48) 



By solving for the first derivatives, a straightforward calculation shows that the 
above set of equations simplifies to 



= — — + m 



dru ' -Sni 7 V 



= + mv + 7 — [vS lf - v k S-t\ , 

dru VK v 



'to 



= 



l s f§ 



\-^ie^f4> S-gS^ SffSg^] , 



(6.49) 

(6.50) 
(6.51) 
(6.52) 



provided that S~s ^ is assumed. Substituting Eqs. HA.30|) . EZt an d then 
Eqs. I|A.34|) . 1A.35() into Eq. (|6.5UI) (see the equations listed in | Appendix A.3\ leads 
to 



o = s A 



- Ck 
Substituting Eq. into Eq. 



JVCr, 



vk c 



Ck v* - v\ 
leads to 



n iv Ck 

d = C 2 = C 3 = C 4 = , C = -J 

7k v k 



1 1 



7 2 2 + j/ 2 



K 



(6.53) 



(6.54) 



implying that S^ — 0, and S*^ , S-^ are constant, from Eq. 1)6. 6|l . and Eqs. I|6.4[l and 
respectively. But from Eq. jOg} it follows that S$ = 0, or A = B = C = D = 0, 
from Eq. ^6.2(1. so that 5^ = and 5-g = as well, from Eqs. (|6.1|l and l|6.3f> 

respectively. This contradicts the assumption S~& ^ so only the case S~x = 
remains to be considered. 
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= 



= 



dru 
dS t 



vS ir - S f4, 



Ca 
— 



dri 



m 



i ^ K r a 
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(6.55) 



(6.56) 



provided that vS if - ^ 0. Substituting Eqs. (|A.30j) . l|6~7ll and then Eqs. (|A.34I) . 



(|A.35|) into Eq. Ij6.56j) we obtain 
= C|(^-4)[44-^]+4co[^%-%]- 7 ^CKC m (^-4)[%-^] .(6.57) 
Substituting Eq. (|6.4|) into Eq. (|6.57|l then gives 

ci = c 2 = c 3 = c 4 = , (6.58) 

implying 



from Eq. 1)6. 6|) . and Eqs. 
satisfied; moreover 







dS, 



rtf) 



dru dru 
and H6.5fl respectively. Hence, Eq. 



(6.59) 
is identically 



(±) C m 71/ (jf 1 2 2\1/ 2 2 \2 1 o 2 r 2r 2 2 yi 

Co ' = —^- — ^[2-u K (l-5i/ K )](u -u K ) +Zv K {v [3 - v K {7 - v K )] 



— + 4(2 + 4) 

7a 



4 n 1/2 



-2- + 4(2 + ^ 2 

7a 



K j 



— 9z/V 



A 



(6.60) 



Next substituting Eqs. (|6.4|l and (|6.5(l and then Eq. (|6.60l) into Eq. I|6.7|l leads to 



= - 



- - ^[i -40. - ^)]|- f 2 4 _ „ 4[1 _ 4(3 + 4)| 

7a 



2z/ 



A' 



+ ^ 2 4[2 - 4(18 - 4)] + 4*4(2 + 4) 



± — {v\\ - 4(3 + v z K )\ + 4(2 + v" K ))\v\\2,v z K + \v z ) - 84] 1/2 \ . (<>.(,! i 

^A 

Finally substituting Eqs. H6.60|l and (|6.61(l into Eqs. I|6.4|) and Ij6.5|l . we obtain 
expressions for the only nonvanishing components of the spin tensor 

m v K v 2 ~ v\ v 2 (2 - Zv 2 K ) + Av\ ± vv K [v 2 {\Zv\ + Av 2 ) - 84] 1 / 2 



4(^-2)-^[i-4(3-4)] 



2 7 Ca 



2 7 Ca v 2 + 2v z K 

m v K {v 2 K ^ - 2) - 4l - 4(3 - 4)]}" 1 



,(6.62) 



v 2 + 2v\ 



{vv K [{v 2 + 2v 2 K ){Av 2 -3-4) 

- 2(^ 2 - 4) 2 ] ± - 34) - 24] [4134 +44 - 84] 1 / 2 } , (6.63) 

which are in agreement with the condition vSf f — ^ assumed above. This 
solution, having constant spin components, was already found in previous work 0. 

Eq. (|6.59|l together with the fact that = 0, S f g — show that the total 
4-momentum P (see Eq. ()3.5Jl ) also lies in the cylinder of the circular orbit 

P = mU + msE, , (6.64) 
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with 



m s = -j—{vS if - UjcSfi) . (6.65) 
It can therefore be written in the form P = fiU p with 

7v r i 1 v + m s /m 7 

= 7 P ef + v p ei , ^ = 7- T^T > ^ = —im + vm s ) , (6.66) 

v 1 + vm s / m 7 p 

where 7 P = (1 — z/ 2 ) -1 / 2 , provided that m + vm s ^= 0. The T supplementary conditions 
can then be written as 

% = . S fi + S f^ v v = . % + = . ( 6 - 67 ) 

the last condition being identically satisfied, and with the equivalent azimuthal velocity 
v p given by 

,(±) - L ^ !%,„,„ + r„2m„ 2 . J- A„ 2 \ - 8,AlV2 



'V = \ 2 To 2 { 3w ^ ± [^(13^ + 4^ 2 ) - 84] V "} , (6-68) 
from Eqs. (|6.62l) and 16.63fl . The reality condition of (|6.68() requires that // takes values 



outside the interval (— v,v), with z> = vkV% V ~~ 13 + 3V33/4 ~ 0.727(/r-; moreover, 
the timelike condition for < 1 is satisfied for all values of v outside the same 
interval. 

From (|6.67l) the spin vector orthogonal to U p is just 7 P _1 S~^Eg. The spin- 
curvature force ()3.7|) turns out to be radially directed 

= 7CI [2% + vS f$ ] e f . (6.69) 
The term on the left hand side of Eq. can be written as 

- — = [mK — m s Ti]ef , (6.70) 
dry 

so that the balance equation 1)3.8)1 reduces to 

m a(U) f = Fj so) + , (6.71) 

where 

ma{U)r = mK , F, (so) = m s n , F- (sc) = 7 (| [2% + uSJ . (6.72) 



7. Conclusions 



Spinning test particles in circular motion around a Schwarzschild black hole have been 
discussed in the framework of the Mathisson-Papapetrou approach supplemented by 
the usual standard conditions. One finds that apart from very special (and indeed 
artificially constrained) orbits where the spin tensor is closely matched to the curvature 
and torsion properties of the world line of the test particle or the static observer spin 
vector is constant and orthogonal to the plane of the orbit, the assumption of circular 
motion is not compatible with these equations. Indeed even in the former case, the test 
particle assumption is then violated except in the limit of massless particles following 
null geodesies, where the spin vector must be aligned with the direction of motion from 
general considerations. The spin-curvature force generically forces the motion away 
from circular orbits, so one needs a much more complicated machinery to attempt to 
study explicit solutions of this problem, solutions which must break the stationary 
axisymmetry. 
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Appendix A. Derivation of the solutions of the equations of motion 

This Appendix derives the solutions of the equations of motion alone Eqs. yi.4| 
and H8.9(l - H3.12(l for the spin tensor along a circular orbit without supplementary 
conditions imposed. Standard elimination and differentiation techniques are used 
to find decoupled second order linear constant coefficient equations for certain spin 
components, from which one may calculate the remaining spin components that do 
not already satisfy decoupled first order such equations. 



Appendix A.l. The v = case 

When the particle is at rest v = relative to the static observers, these equations 
reduce to 

0=-^-"kCkS { }, (A.l) 
dSxi 

(A.2) 

= — ^ - vkCkSm , (A.3) 

<±Tu 

= -t—+vk^k-t Jl , (A.4) 
dru dru 

= - 2 & S * + mVK ^K , (A.5) 



'U 

d 2 ^ , dS., 

° < A - 7 > 

Eq. (|A.2|) implies that Sg^ = c\ is constant. Solving Eq. (|A.3|I for dS f §/dTjj, and 
substituting the result into Eq. (| A. 6|> leads to the decoupled equation 



_ d2s to , w _ _ M(r-3M) 

Similarly solving Eq. (|A.1|) for dS^/dry, and substituting the result into Eq. (|A.7|I 
leads to an analogous decoupled equation 

= ^?+-o^- (A.9) 

Eq. (|A.4f) leads immediately to m = —i^kCkS^ + c m and substituting this into (|A.5|I 
yields 



0= ^+<4S tf + VK CKC m , = Ck(2 + 4) 1/2 = J ^ 3 ( ( 2 ;_ ^ .(A.10) 

The three second order constant coefficient Eqs. (|A.8ll . (| A. 9|> and (|A. 10|l are easily 
integrated, from which expressions for the remaining components of the spin tensor are 
then obtained from Eqs. (|A.1(1 and i|A.3|) . These have either oscillatory or exponential 
solutions depending on whether the squared frequencies in Eqs. HA.8(1 and (|A.10(1 are 
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positive or negative, or linear solutions when zero. This distinguishes the two intervals 
2M < r < 3M and r > 3M, whose corresponding solutions are given by Eqs. jSJ 
and l|4.3[) respectively. The special case r = 3M can be easily discussed by setting 
vk = 1 and (k — \^3/(9M) in Eqs. HA.1|) - (|A.7|I . The corresponding solution is given 
by Eq. g^J. 



Appendix A. 2. The v = ±i>k case 

When the particle moves along a geodesic with v = ±vk, Eqs. I|3.2JI - [I3.4JI and (|3.9|1 - 
(|3.12[) simplify to 

dS xl 



dru 
dS §4> 



T vk 



ds if 

dru 

dS{§ (k 



dT V 7 X 



dS, 



;■() 



dru 
d 2 S. 



±UlK[S S iTV K S iS ] 



tei > 



1 dm 
±— — ±Kk1k 



dr^ v K dru 



do ^ 



dS*i 



dr v 



T vk 



dru 



d 2 S. 



dr, 



— T v k- 



d2 % n CK*Si 



t</> 



dr?j jk dru 



T v K ^ + &[S a ± , 



= 



d 2 S. 



Lip 



dm 
dru 



± fit- — ± 2(^7^ 



d% 
dTy 



dS. 



T v K 



dru 



(A.ll) 
(A.12) 
(A.13) 
(A.14) 
(A.15) 
(A.16) 
(A.17) 



The difference of Eqs. I|A.14|) and (|A. 17|) leads to dm /dry = 0, so m — c m is 
constant, which then implies from the same equations that 







d 2 S. 



10 



dr 2 



± 2( K IK 



dS {f 



dS £l 



dm 



T vk 



dm 



Integration of Eq. (jA.llfl yields 

% = ±v K S if + ci , 
and using this to replace S-^ in Eq. (|A.18(1 and then integrating gives 



(A.18) 
(A.19) 

(A.20) 



dry 7^ 

Using Eqs. (|A.19I) and HA.20I) to replace and dS^/diy in Eq. (|A.15(I then leads 
to the decoupled second order equation 

= ^f+^S if ±2 lK (KC 2 , c 2 = 0K4-37I) 1 / 2 = ^ ( ( ;: 3 6 g .(A.21) 

Taking the tj/ derivative of Eq. (|A.12|) and using it and Eq. (|A.13j) in Eq. I|A.16(1 leads 
to the second order equation 

d 2 S tl 







dr 2 



u 3 S t§ + 3(k1kVkS ss> , uj 3 = C 



K 



M 



1/2 



(A.22) 
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Finally using Eq. I|A.22[1 and Eq. (|A.13|I in the equation obtained by taking the tjj 
derivative of Eq. (|A.12|) gives a second decoupled second order equation 







W4 



2) 



1/2 



M 



r 3 (r-3M) 



(A.23) 



Integrating the two decoupled second order equations l|A.21() and (jA.23fl . one can 
then integrate Eq. (|A.22|) too. The remaining components of the spin tensor are then 
determined by Eqs. (|A.13|) . i|A.19(l . and i|A.20() . Note that the frequencies lo^ and 0J4 
agree only for jk = 1, or Vk = 0, which would imply M = 0: so this special case is 
not relevant. 

Now the two intervals 3M < r < 6M and r > 6M have differing signs for the 
squared angular velocities, and the corresponding solutions are given by Eqs. I|5.1[) 
and H5.3J) respectively. The special case r = 6M can be easily handled as well, setting 
v K = 1/2 and Ck = \/6/(36M) in Eqs. I|A.11HA.17|) . Tne corresponding solution is 
given by Eq. (|5.2jl . 



Appendix A. 3. The general case: v 7^ and v 7^ ±vk 

Here we deal with the general form of Eqs. (|3.2() - (|3.4() and (|3.9|) - (|3.12() . Solving 
Eqs. I|3.2|) - (|3.4|) for their first terms and substituting these derivative terms into 

Eqs. EUt-dlllll, one finds 

d 2 S 12 



t<t> 



dr 2 

d 2 Sjf 



v vk aru 



3 frr 2 
71/ IV 

VK 



dr f , 



^ v d Tu 
■ - T [(3v K -iy-2v K ]St f 



.A Cff 1 



K 



- S K )S fi + m 7 ^(v 2 v 2 K ) , 



VK 



d 2 S t 



7 3 C 



^ ,,2 



cLS' 



;"0 



+ 7 4 *#[(1 + 2^ 2 - 3^]5^ 
Ik v k 



dr 2 v 



d 2 S; 



tcj> 



dr 2 v 



Ik v k 

CK r 2 
: 7^ ll/ 

VK 



dr r , 



2] d% 
drr 



1 2 ^{v 2 -Vk)[v 2 + v k -1]S u 

V T.' 



dm 
dr r . 



C7 Vk t utc/ 

Solving Eqs. (|A.24|I and (|A.27|1 for dm/drjy and d 2 S^ i/dr 2 ,, one finds 



dm Ck , 2 2 \ 



dry 

*¥ 

dr 2 

U7 [/ 



dry 



7^— -Sf, 



= - 2 



'if /,,2 



7 v (k dS if 7 2 Cj 



Solving Eq. (f^T^|) for S^g, leads to 

Q _ 1 V K 1 



J t<j> 



7 Cif v 2 - v 2 K 



dSt. 



dS £ - 



dru dru 



Using Eq. (|A~30|) in Eq. (|A~28)l leads to 



Ck r „ 



A' 



(A.24) 

(A.25) 
(A.26) 

.(A.27) 

(A.28) 
(A.29) 

(A.30) 
(A.31) 
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,d 2 S. 



d 2 S. 



7 2 C 



= [(1 - 2viy + ^ K ]^f + v[v 2 + A ~ 2]-=^ + W ~ -k) 2 S, 



Ik v K 

vCk r,,2 ,.2 \|71 i„2\„2 i ,2 In , , ,2 \i c i ^ A '/,.2 



+ r W - «4)[(1 - 4*4K + *4(2 + ^)]% + - ^)'c m . (A.32) 



A" 



^A 



Using Eq. (fO0|) in Eq. (fQ9l leads to 

d2 % „d 2 % T^Ci 
7a 4 

Then solving these last two equations for d 2 ^. i/dr^ and d 2 £^/dTy gives 



(A.33) 



d 2 % 



^Ca 
7 4 



+ 7V 



4(2 + ^) 



f 0' 

-y 2 r 2 

7a u K 



3 Ca / 2 2 \ 
2 ,,2 " C + 7 — - ^a) c « 



'A 



K 



(A.34) 



d 2 5. 



dr 2 



4J(^ 2 + + 7 2 K^(1 + 2*4)% 

7a ^a 

2 z/ 2 (l-2^ 2 ,) + i/| f s„Ci?,2 



c 4- Tv — [v - v K )c m , 



whose solution is given by Eqs 
Eq. HA.30|1 . one obtains l(6~l)|) . 

Substituting Eqs. I|3.3|l and (|3.4|l into Eq. (|3.11[) gives 

d 2 &, 



(A.35) 

and Q6.5J1 . Then substituting these solutions into 







'10 



dT 2 |^ + 37 2 K|^, 



7a 



(A.36) 



Taking the tjj derivative of Eq. I|3.3|l and using Eqs. (|A.36I) and (|3.4I) . one gets 



d 2 S 1 -- 
dr 2 ^ 



- ->"7v 

^A 



(A.37) 



This is easily integrated to give the solution 1(6.1(1 . The remaining components of the 
spin can be obtained directly from Eqs. I|A.36|) and l|3.4|l yielding Eqs. 1(6.2(1 and (|6.3|) . 

The frequencies Q and f^i agree for the particular value of the azimuthal velocity given 

by 



u = ±^(l+2u 2 K )- 1 / 2 
7a 



(A.38) 



The solution corresponding to this special case is given by Eqs. 116.12(1 - 1)6.14(1 together 
with Eqs. ((6.4(1 - 1(6.7(1 evaluated at v = vq. 



References 

[1] Mathisson M 1937 Acta Phys. Polonica 6 167 

[2] Papapetrou A 1951 Proc. Roy. Soc. London 209 248 

[3] Corinaldesi E and Papapetrou A 1951 Proc. Roy. Soc. London 209 259 

[4] Pirani F 1956 Acta Phys. Polon. 15 389 

[5] Tulczyjew W 1959 Acta Phys. Polon. 18 393 

[6] Bini D, de Felice F and Geralico A 2004 Class. Quantum Grav. 21 5427 

[7] Bini D, Carini P and Jantzen R T 1997 Int. J. Mod. Phys. D 6 1 



Spin precession in the Schwarzschild spacetime: circular orbits 

[8] Bini D, Carini P and Jantzcn R T 1997 Int. J. Mod. Phys. D 6 143 

[9] Iyer B R and Vishveshwara C V 1993 Phys. Rev. D 48 5721 

[10] M0ller C 1949 Commun. Dublin Inst. Adv. Studies A 5 3 

[11] Mashhoon B 1975 Ann. Phys. 89 254 

[12] Suzuki S and Macda K 1997 Phys. Rev. D 55 4848 

[13] Semerak O 1999 MNRAS 308 863 

[14] Hartl M D 2003 Phys. Rev. D 67 024005 

[15] Hartl M D 2003 Phys. Rev. D 67 104023 



